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Abstract 

We report the experimental realization of the purification protocol for single 
qubits sent through a depolarization channel. The qubits are associated with 
polarization encoded photon particles and the protocol is achieved by means 
of passive linear optical elements. The present approach may represent a 
convenient alternative to the distillation and error correction protocols of 
quantum information. 
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Modern quantum data processing using realistic (imperfect) quantum gates and long- 
distance quantum communication in the presence of a noisy environment requires a large 
supply of qubits with a high degree of purity. Indeed the fidelity of most quantum informa- 
tion (QI) protocols critically depends on the preservation of the purity of the QI carriers. It 
is therefore crucial to develop techniques that protect quantum states from the unavoidable 
losses and decoherence processes accompanying the transmission. One of these techniques is 
the quantum error correction [1,2] which works by encoding the quantum state into a higher- 
dimensional Hilbert space. Alternatively, one can distribute several copies of an entangled 
state and extract fewer highly-entangled states by means of entanglement distillation [3-7] 
in order to be able to subsequently transmit an arbitrary state with high fidelity by quantum 
teleportation [8-11]. Yet another option is to transmit several copies of the state over the 
noisy channel and then purify the resulting mixed states at the receiver's station [12-14]. 

The present work realizes the purification procedure that was theoretically proposed by 
Cirac et al. in 1999 [12]. It addresses the issue of the purification of iV equally prepared 
qubits in the mixed state p — £ \ 4>) (0| + |(1 — where < £ < 1. This procedure allows 
to distill from a set of mixed states a subset of states with a higher degree of purity, i.e. it 
probabilistically increases the purity by filtering out some of the noise. The procedure is 
based on a set of projections onto the symmetric subspace of the N qubits (i.e. the subspace 
spanned by all the states that are invariant under any permutation of the N qubits) and onto 
orthogonal subspaces that contain symmetric subspaces for subsets of the initial N qubits. 
The procedure is designed to be optimal and universal, i.e., it acts with the same fidelity for 
all input states. Since it is optimal, the purity cannot be further increased by any means. 
In this paper we consider the case of two qubits, i.e. N = 2. The purification procedure for 
N = 2 reduces to a projection of the two-qubit state onto the symmetric subspace, and it 
is equivalent to the symmetrization procedure proposed as a theoretical method to stabilize 
quantum computation in the presence of noise [15]. 

For N = 2 the purification procedure works as follows. Consider two independent qubits, 
a and b, both originally in the state that are transmitted over a noisy channel from which 
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they emerge in a mixed state represented by the density matrix p a = £ \<p) (<p\ + |(1 — £)I = 
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where is a state orthogonal to Our goal is then to 
purify the transmitted qubits in order to obtain two qubits that are as close as possible to 
the original state \4>). The overall 2-qubit input state p™ b = p„™ <g> p l b n is expressed in the 
basis {\<f>) a \4>) b M) a \t^) b , |^-> a \4) b ,\^) a |0 ± >J by the matrix 
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As mentioned above, the purification protocol consists of the projection of the 2-qubit state 
onto the symmetric subspace: if the projection is successful we obtain two equal output 
qubits that are the optimal "purified" ones, otherwise we discard the output states. We 
note that this protocol can be implemented, for every qubit encoding, by a quantum circuit 
requiring an ancilla qubit and a Toffoli gate [16]. After a successful projection the output 
qubits are in the state 
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where II = I & — \^ab)(^ab\ * s the projector onto the symmetric subspace and \^ab) = 
(|01) — |10))/\/2 is the singlet state of two qubits. The success probability of the procedure 
is p = TrfUp^lT] = ^p-. Since p™£ belongs to the symmetric subspace, the reduced density 
matrices of the resulting single qubits, expressed in the basis {]</>) , are found to be 

identical, 
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where £p = £ > £ and the purification gain factor is 77 = 3^2 ■ Note that p and i] are 
related by the equation r\p = 1 so a higher purification gain factor is necessarily accompanied 
by a lower probability of success. 

We report the implementation of the above protocol for qubits encoded in the polarization 
of single photons (see Figure 1). The qubit to be purified is ^ \4>) {4>\ + ^ [(j) 1 ) where 
|0) = a \H) +/3 \V) and \H), \V) respectively correspond to the horizontal and vertical linear 
polarizations. In the present experiment, pairs of photons with wavelength A = 532 nm and 
coherence time r co h = 80 fs, were generated in a Type I, BBO crystal slab in the product 
state \H) a \H) b by spontaneous parametric down conversion (SPDC) process excited by a 
CW fourth-harmonic-generation laser (Coherent Verdi +MBD-266). The output state was 
first encoded in the state \4>) a \4>) b by means of two equal waveplates (wp) WP(\4>)) and 
then each photon, injected into a noisy channel P, emerged in the mixed state: p a = p b 
= £ 10) (01 + (1 — £)§■ The mixed qubits, associated with the two modes k a and k b , were 
linearly superimposed at beam-splitter BS with a mutual time delay At micrometrically 
adjustable by a translation stage with position settings Z = 2 Ate, with c denoting the 
velocity of light. The value Z = was assumed to correspond to the full overlapping of 
the photon pulses injected into BS, i.e. to the maximum photon interference leading to 
the simultaneous detection of two photons on either output modes k-y or k 2 of BS [17] 
[18]. Recently it has been shown that the projection of the overall state in the symmetric 
subspace, precisely the one implying the present purification procedure, is unambiguously 
identified by the maximum interference condition: Z = [18] [19]. 

Let us give more details about the realization of the two equal depolarizing channels P 
and P' operating on the BS input modes k a , k b , respectively. Each channel implemented the 
quantum map S (p) = £p + (1 — ^)£dep (p) where £dep (p) maps any unknown input state 
p into a fully mixed one. This transformation can be achieved by stochastically applying the 
full set of Pauli operators {I, a x , a y , a z } with the same statistical weight, that is, Edep (p) = 
\(IpI + o~ x po~x + °~yP°~y + o~zPO~z)- Let us consider here only one of the noisy channels, say 
P, the one operating on the mode k a . The £ (p) map was realized by means of a pair of 
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equal electro-optic (EO) LiNb0 3 Pockels cells (P-cell), Px and Pz, carefully aligned with 
a 45° mutual spatial orientation of the optical axes of the EO crystals (see Fig. 1). All 
P-cells were Shanghai Institute of Ceramics devices with a ^-voltage: Vx = 390 V. Each 

A 2 

P-cell of the pair was driven by a CW periodic square-wave electric field with maximum 
V = Vx, fixed frequency / = T _1 and variable pulse duration r corresponding to a duty- 
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cycle v = t/T adjustable in the range: < v < 1/2. The excitation pulses feeding the two 
P-cells were mutually delayed by a time equal to r/2 (see inset of Fig. 1). Consider a single 
excitation cycle. In the time intervals At = r/2 in which only one P-cell was active, either 
the a x or the a z transformation was implemented depending on the corresponding crystal 
orientation. In the interval At = r/2 in which both P-cells were simultaneously active, the 
o y transformation was realized. In summary, each operators I, a x , a y , a z was applied to the 
input state over a time At = r/2 and the total depolarizing process lasted a time 2r over 
each period T, thus achieving an average depolarizing fraction (1 — = 2u. In order to avoid 
any correlation between the two qubits to be purified, p a and p b , the two channels P and P' 
were driven by different frequencies: / =6 KHz and /' = 1.7 x /. Correspondingly, different 
values of r were adopted for the two channels in order to realize, within each experimental 
run, equal values of £ for the two input qubits. 

In the analysis we have assumed an identical preparation of the two input qubits, while 
the output ones are described by the same density matrix p = p a = p b . With this assumption, 
carefully checked over each channel, the verification of the purification procedure lies on the 
tomography of the density matrix of one of the input and one of the output qubits. For the 
sake of simplicity, we only analyzed the measurements performed on the BS output mode 
ki (see Fig. 1), selecting counts in coincidence between the detectors [Di,D 2 ] to trigger the 
realization of the projection of p onto the symmetric sub-space. The detectors Di,D 2 were 
coupled to mode k-y by a 50:50 beam-splitter BS\. D\ provided the measured outcomes of a 
simple tomographic setup consisting of a A/2-wp, a A/4-wp and a polarizing beam splitter 
(PBS). 

Consider first the projector switched off, by setting Z ^> CT co h, i.e., by spoiling any 
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interference on the photons impinging on BS. A tomographic reconstruction of the qubit in 
the mode k[ based on the measurement of the corresponding Stokes parameters by 4 different 
settings of the wp's |, j was undertaken. It is easy to see that this qubit, corresponding to 
the qubit to be purified, is expressed by the density matrix p 1 = - 2 b [20]. By turning 
on the projector, i.e., by restoring the BS interference setting Z = 0, the mode k\ contains 
the two photons described by the density matrix p°^ . In this case, we measured on the 
mode k[ the purified qubit p\ = p° ut . From the density matrices reconstructed in absence 
and in presence of interference, we obtain £, £ p and thus the purification factor r\ = £ p /£. 
In addition, from the coincidence rates determined for Z = and for Z cr co h we inferred 
the success probability p of the purification protocol. We may check that an increase of the 
purification gain factor for any qubit pair, i.e., a larger rj, corresponds to a lower success 
probability of the overall protocol, as expected. In Fig. 2, we plotted the experimental 
values of rj and p obtained for different £'s, i.e., for different experimental values of v — 
(1 - 0/2, for three input states: \H), \L) = 2~*(\H) + \V)), and \E) = [cos(37r/16) \H) + 
isin(37r/16) \V)] corresponding, respectively, to horizontal, 45°-diagonal, and a very general 
elliptical polarizations of the input qubits. The mutual agreement of the data for different 
input states demonstrates the universality of the purification procedure. The deviations of 
the experimental data from the theoretical values were mainly due to the imperfections of 
the optical components. In particular, the non-ideal properties of the main BS were found 
to be highly critical. In order to achieve the projection onto the symmetric subspace, the BS 
transmittances T# and Ty for the H and V polarization modes should be equal, with a high 
level of precision, and any difference between T H and T v partially spoils the purification. 
Notice, however, that deviations of T H = T v from 50% only decrease the success probability 
p but do not alter the purification gain factor 77. 

We may generalize the above method by accounting for any possible asymmetry of 
the preparation of the input qubits. Allowing the input qubits to have different degree 
of mixedness, i.e. = C 10) (01 + (1 — C)f> PT — K \4>) (01 + (1 — the output 

qubits are still in the state given by Eq.3 with £p=2(£ + re)/(3 + kQ. The purification 
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factor is 77 = £p/£ = = 1/p where £ = |(C + K ) is the average input mixed- 

ness factor. This process may be investigated by recourse to the quantum "relative en- 
tropy" that measures the closeness of any output state a with respect to a correspond- 
ing input pure state p, e.g. after propagation through a noisy communication channel: 
S(p || a) = Tr(pLogp) — Tr(pLoga) [21]. Suppose that two qubits are equally prepared in 
the pure state p = \4>) a {<f>\ a <g> \<p) b (<p\ b , S(p) = . After corruption by noise the entropy is: 
S(p || p rn )=Log [|(1 + £)] +Log [|(1 + re)] . If the qubits are further purified by symmetriza- 
tion the following result is obtained: S(p \\ p out )=Log [|(1 + £)] + Log [|(1 + re)] — Log?]. 
Then the symmetrization leads to the positive information gain AS=S(p || p m ) — S(p \\ p out ) 
= Logrj at the expense of a reduced rate p of success: AS = —Logp, | < p < 1. 

An interesting case is represented by the purification of a mixed state by a jrare state, 
e.g. by the initial conditions (=1 and re=0. This precisely corresponds to the probabilistic 
quantum cloning process recently realized by our Laboratory in Rome by a symmetrization 
procedure [16], [18]. There £ = |, and a purification gain factor 77 = 4/3 was attained with 
a success probability: p = 3/4. 

In conclusion, we have experimentally demonstrated the optimal purification of two 
depolarized qubits using the interference of two photons at a beam splitter, conditionally 
effecting symmetrization. The experimentally observed purification gain factors are in very 
good agreement with the theoretical estimates. We therefore envision that single qubit 
purification may become a viable procedure for protecting quantum states against noise. 
Note also that the projection on the symmetric subspace of more than two qubits can 
be carried out with a sequence of beam splitters, which is in the reach of present optical 
technology. 
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Figure Captions 

Figure 1. Experimental setup for the optical implementation of the purification procedure. 
Inset shows the realization of the depolarizing channel employing two Pockels cells. 

Figure 2. Experimental results of the purification procedure for different input qubits 
corresponding to the encoded polarizations: [if), \L) = 2" 1//2 (\H) + \V)), and \E) = 
(cos (|) \H) -Msin (|) \ V)) with 9 = |7r. Filled markers denote the experimental purifi- 
cation factor rj data while open markers denote the experimental data of the procedure 
probability p. For simplicity we report only one error-bar for the probability values ob- 
tained. The statical error is the same for all data reported. 
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